Submitted to ApJ 

Weak Lensing Measurements: A Revisited Method 
and Application to HST Images 

Jason Rhodes 

Physics Department, Princeton University, Jadwin Hall, P.O. Box 708, Princeton, NJ 08544; 

jrhodes@pupgg.princeton.edu 

and 

Alexandre Refregieif] 
Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544; 

r efreg@ astro . pr inceton . edu 

and 

Edward J. Groth 

Physics Department, Princeton University, Jadwin Hall, P.O. Box 708, Princeton, NJ 08544; 

groth@pupgg.princeton.edu 

ABSTRACT 

The weak distortions produced by gravitational lensing in the images of background 
galaxies provide a unique method to measure directly the distribution of mass in the 
universe. However, because the induced distortions are only of a few percent, this 
technique requires high precision measurements of the lensing shear and cautious 
corrections for systematic effects. Kaiser, Squires, & Broadhurst (1995) proposed 
a method to calibrate the ellipticity-shear relation in the presence of Point Spread 
Function (PSF) anisotropics and camera distortions. Here, we revisit the KSB method 
in the context of the demanding search for weak lensing by large-scale structure. 
We show that both the PSF and the camera distortions can be corrected for using 
source moments, as opposed to ellipticities. We clarify the applicability of some of 
the approximations made in this method. We derive expressions for the corrections 
which only involve the galaxy moments. By decomposing the moments into spinors, 
we derive an explicit relation between the shear and the average ellipticity. We discuss 
the shortcomings of the method, and test its validity using numerical simulations. As 
an application of the method, we repeat the analysis of the HST-WFPC2 camera 
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performed by Hoekstra et al. (1998). We confirm the presence of sizable (~ 10%) 
PSF ellipticities at the edge of the WFPC2 chips. However, we find that the camera 
distortion is radial, rather than tangential. We also show that the PSF ellipticity 
varies by as much as 2% over time. We use these measurements to correct the shape of 
galaxies in the HST Survey Strip ("Groth" Strip). By considering the dependence of 
the ellipticities on object size, we show that, after corrections, the residual systematic 
uncertainty for galaxies with radii greater than 0.15 arcsec, is about 0.4%, when 
averaged over each chip. We discuss how these results provide good prospects for 
measuring weak lensing by large-scale structure with deep HST surveys. 

Subject headings: cosmology: observations - gravitational lensing - methods: data 
analysis - techniques: image processing, photometric 

1. Introduction 

Weak gravitational lensing produces coherent distortions in the images of background 
galaxies. This effect provides a unique method to measure directly the distribution of mass in the 
universe (for reviews see Schneider et al. 1992; Narayan & Bartelmann 1996; Mellier 1999). This 
technique is now routinely used to map the mass of clusters of galaxies (for a review, see Fort & 
Mellier 1994). A search for weak lensing by large-scale structure is the subject of much recent and 
on-going theoretical and observational effort (eg. Villumsen 1995; Stebbins, et al. 1995; Kaiser 
1996; Schneider et al. 1997; Van Waerbeke et al. 1998; Refregier et al. 1998; see Refregier 1999 for 
a bibliography). The main difficulty lies in the fact that the lensing distortions are small (~ 10% 
for clusters and ~ 1% for large-scale structure), thus requiring high precision measurements and 
tight control of systematic effects. 

Kaiser, Squires and Broadhurst (1995, KSB) have developed a method to correct for the 
major systematic effects, namely the anisotropy of the Point-Spread Function (PSF) and camera 
distortions, and to calibrate the relation between galaxy ellipticities and lensing shear (for other 
methods, see also Bonnet & Mellier 1995; Schneider & Seitz 1995). Further elements of their 
method were presented in Luppino & Kaiser (1997, LK) and in Hoekstra et al. (1998, HFKS). 

Recently, Kaiser (1999) pointed out that the KSB method had several shortcomings, all 
stemming from the fact that most PSFs encountered in practice are not sufficiently compact. 
Kaiser then proposed another method based on the explicit construction of the post-convolution 
shear operator. Another alternative method was recently proposed by Kuijken (1999). In this 
different approach, a sheared and convolved isotropic model is fitted to the galaxy image, so as to 
derive an estimator for the shear. These two methods are promising, but both require complete 
knowledge of the 2-dimensional PSF function, while PSF measurements are sparse due to the 
finite number of stellar images. 
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Here, we revisit the KSB method, which has the advantage both of being linear and of relying 
only on the first multipole moments of the PSF and galaxy images. We focus on the demanding 
search for weak lensing by large-scale structure. We show that both the PSF and the camera 
distortions can be corrected for using source moments, as opposed to ellipticities. We clarify 
the applicability of some of the approximations made in this method, and show how the weight 
function for stars can be chosen to be different from that for galaxies. We derive expressions for 
the corrections in term of the moments only. By decomposing the moments into spinors, we derive 
an explicit relation between the shear and the average ellipticity. We discuss the shortcomings of 
the method discussed in Kaiser (1999) and Kuijken (1999), and test its validity using numerical 
simulations. 

As an application, we consider weak lensing measurements with the Hubble Space Telescope 
(HST). The small PSF and absence of atmospheric seeing makes HST an ideal instrument for weak 
lensing measurements (eg. Kneib et al. 1996; HFKS). We reproduce the analysis of HFKS who 
studied, in detail, the PSF and camera distortion of the WFPC2 camera onboard HST. We apply 
these calibrations to the galaxies in the HST Survey Strip known as the "Groth Strip" (Groth et 
al. 1995, Rhodes 1999). In particular, we show how the PSF anisotropy and camera distortions 
depend on the size of the galaxies. We also discuss the prospects of weak lensing measurements 
with HST. A description of our search for weak lensing by large-scale structure with the Survey 
Strip will be presented in Rhodes, Refregier, &, Groth 1999 (see also Rhodes 1999). 

Because the method is somewhat complex, we provide a practical summary in §[2] in which 
we point to the results and equations which are of direct practical interest. In §||, we describe 
measures of object shapes. Next, in §|] and §|BL we show how these measures are affected by 
the two main classes of deformations, namely convolutions and distortions, and derive explicit 
expressions to correct for them. In §[| we study the effect of shear combined with a convolution 
on the observed ellipticity. In §^, we discuss the shortcomings of the method, and test its validity 
using numerical simulations. In §||, we apply our method to HST observations. In particular, we 
consider measurements of the camera distortion, globular cluster observations, and finally apply 
our results to the Survey Strip. Our conclusions are summarized in §||. 

2. Overview of the Method 

The purpose of this method is to provide a measure of the shapes of galaxies, to correct for 
instrumental effects, and to derive an estimate for the weak lensing shear. The images of galaxies 
are assumed to be altered by two kinds of operations: distortions (sometimes simply called shear), 
which are intensity-conserving mappings between the source plane and the image plane (see 
Eq.Jl^|), and convolutions (or smear; see Eq. j36|). Specifically, we follow KSB and HFKS and 
assume that a galaxy image is altered by a distortion due to gravitational lensing, followed by a 
convolution by a generally anisotropic PSF, and then by a distortion due to the camera optics. 
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Schematically, we can write 

intrinsic image — > lensing(<9) — > PSF(*) — > camera(<9) — ► observed image, (1) 

where d and * symbolize distortions and convolutions, respectively. Here, we revisit the KSB 
method. We show how to correct for the two latter effects and how to calibrate the sought-after 
effect of lensing. 

The synopsis of our method is shown on figure |]. Since the method is somewhat complex, we 
provide here a summary of the method, and point to the equations which are of direct practical 
interest. References to these equations can also be found on this figure. 

The basis of the method is to characterize the shapes of galaxies and stars by measuring their 
multipole moments «7y, Jijki, etc. (Eq. ||). To enforce convergence in the presence of image noise, 
these moments are measured with a weight function w(9) which we choose to be a gaussian of 
width lo (Eq. ||). This choice allows us to write all the expressions for the corrections in terms of 
the moments. Unlike KSB, we perform all the corrections using moments, and postpone the use 
of ellipticities (Eq. p2| ) until the last step. This has the advantage of keeping track of the size of 
the PSF and of the galaxies, and thus of reducing the noise resulting from the corrections, if the 
PSF size varies across the field. 



The first step consists of deriving the distortion matrix ( ^?? mera (Eq. IfLql ) for the camera 
distortion. HFKS showed that this could be achieved by considering astrometric shift solutions, 
such as that of Holtzmann et al. (1995) for the WFPC2 camera on HST. The resulting shear 



pattern for WFPC2 is shown on figure y. (See discussion in §8^. about the difference between this 
figure and the results of HFKS). 

In the second step, we derive the PSF moments from stellar images. For this purpose, the 
stellar multipole moments J*j and J*j ki are measured using an optimally chosen weight-function 
width w*. These moments are then corrected for the weight function to provide (an approximation 
to) the unweighted moments of the PSF (Eqs. [p^j55| ). This allows us to use a different 
weight-function width to for the galaxies, and thus to improve the sensitivity. The PSF moments 
are then corrected for the camera distortion (Eqs p4| , p5"| with uj — > oo). This provides the 
corrected unweighted PSF moments Pjj and Pijki (Eqs. (3^j5l|), that can be then be interpolated 
across the field using low-order polynomial fits for each component. Figures [8| and || show the 
resulting ellipticities of the WFPC2 PSF derived from two globular clusters, while figure 10 shows 
the interpolated PSF ellipticities derived from a combination of the two after the moments have 
been corrected for weighting and camera distortion. 

In the third step, we measure the galaxy moments and correct them for instrumental effects. 
The galaxy moments Jj,- and Jijti are measured using an optimal weight-function width uj. 
They are then corrected for the camera distortion using ^>9|' mera (Eqs. [ p^j25[| ). The PSF can be 
decomposed into an anisotropic and an isotropic part (Eq. p3|] ). We correct the galaxy moments 
Jij for the PSF anisotropy, by expanding in powers of the PSF ellipticity (Eq. |46|] ). Since, by 
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construction, the isotropic part of the PSF is a gaussian (Eq. p4| ), we can perform the isotropic 
correction exactly for Jij (Eq. p9| ). The fourth order moments Jijki are then approximately 
corrected for the PSF (Eq. |5(|). This provides us with the corrected galaxy moments J^ 9 and 
Jijki which are now effectively weighted by the new weight-function width ui g (Eq. [p/^]). 

The final step consists of measuring the weak lensing shear by averaging over an ensemble 
of galaxies in a region of the sky. For this purpose, we compute the ellipticity q of each galaxy 
from its corrected moments J^ B (Eq. fT2|). The weak lensing shear ji is then computed from 
the average ellipticity (e^) (Eq. |3(|). In this last step, one should remember to use uj g rather 
than u as the weight-function width. The relationship between 7j and (e^) is greatly simplified by 
considering the rotation properties of the multipole moments (see § |3.2|) . Note that our method 
avoids the complications of the pre-smear and post-smear shear susceptibilities discussed in LK 
and HFKS. In §|8[ we discuss the application of our method to the HST Survey Strip. 



3. Source Shape Characterization 

In this section, we show how object shapes can be characterised using multipole moments and 
related quantities. We also study the rotational properties of the moments and decompose them 
into spinor representations. 



3.1. Moments 

Let us consider a source with intensity i{9). As a first step, we find the centroid 9° of the 
source by solving 

' d 2 9 {0i - 9?)w{9 - 9°)i(9) = 0, (2) 



where w(9) is a weight function introduced to ensure convergence in the presence of noise. In this 
paper, we will consider a normalized gaussian weight function, 

»(«)^^, (3) 

which has convenient analytical properties. In practice, equation (|2|) can be solved iteratively, 
by fixing the weight-function width to to an initial estimate of the source size (see § |8,3| for a 
description of our choice of uj for the HST survey strip). 

To characterize the source shape, we consider the weighted multipole moments of the source 
intensity, 

I = J d 2 9 w{9)i(9), = J d 2 9 9i9 jW {9)i(9), I ijk = J d 2 9 9 i 9 j 9 k w{9)i(9), etc. (4) 

where the origin of the coordinate system was chosen to coincide with 9°. (By definition, the 
dipole moment ij vanishes in this coordinate system). It is usually more convenient to consider 
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the normalized moments 

Jij — Jijk = Iijk/I} etc.. 

The normalized quadrupole moments J™ can be diagonalised as 

' a 2 ^ 



J = R (-a) 



b 2 



R(-a), 



(5) 



(6) 



where a and b are the (weighted) major and minor radii, a is the position angle measured 
counter-clockwise from the positive x-axis, and T stands for the transpose operation. The rotation 
matrix R is defined as 

cos(y) — sin(</?) \ 
sin((/p) cos(ip) I 



(7) 



Inverting this relation yields 

{a\b 2 ) 
tan 2a 



Jll + J22±V( Jl1 - J22 ) 2 + 4J 12 
2Jl2 



J 11 — J22 



(8) 
(9) 



3.2. Rotational Properties 

To study the rotational properties of the tensors defined above, let us consider a new 
coordinate system which is rotated counter-clockwise by an angle ip from the original positive 
x-axis. In this new coordinate system, the components J-j of the normalized quadrupole moments 
are related to the unrotated components Jij by J[a = Rik(—^p)Rji(—^p)Jki, and similarly for tensors 
of higher order, where R is the rotation matrix defined in equation (^). In general, any tensor in 2 
dimensions can be decomposed into 2-component (or 1-component for scalars) £-spinors u, which 
rotate as 

v'^Riji-l^Vj, £ = 0,±1,±,2,..., (10) 

under this change of coordinates. (These spinors form irreducible representations of £0(2), the 
rotation group in 2 dimensions). 

For instance, the quadrupole moment Jij, which, being symmetric, consists of 3 independent 
components, can be decomposed into a scalar 

d 2 ^\(J ll + J 22 )= l -(a 2 + b 2 \ (11) 

which is the mean square radius of the source, and a spin-2 tensor 

_ { Ju — J 2 2, 2J 12 } a 2 - b 2 , . . 

e l = = {cos(2a ,sm (2a)} (12) 

Ju + J22 a z + b A 
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a spin-2 tensor, 



which is the ellipticity and has been normalized to follow the weak-lensing nomenclature. The 
component t\ (€2) corresponds to stretches and compressions parallel to (at 45° from) the 
coordinate axes. 

By considering infinitesimal rotations, we can also decompose the fourth-order moment Jijf-i, 
which has 5 independent components. We find that it can be decomposed into a scalar, 

A = (J mi + 2J 1122 + J 2222 ) /(2dV), (13) 
Ml = (-Jnn + J2222)/{2d 2 Lo 2 ), 

M2 = -2(Jni2 + Ji 222 )/(2d 2 ^ 2 ), (14) 
v\ = (Jim ~ 6Jn 22 + J2222)/(2d 2 ^ 2 ), 

V2 = 4(Jni2 - Jl222)/(2^ 2 ). (15) 

For future convenience, these spinors have been normalized with the weight-function width u 
(Eq. H) and the scalar d. These decompositions are useful for simplifying tensors which are 
averaged over an ensemble of randomly-oriented galaxies (see §4.4). 



and a spin-4 tensor, 



4. Distortion 

In this section, we study the effect of distortions on the source moments. In practice, 
distortions arise from the instrument optics and from weak lensing. We show how the former can 
be corrected, and how the latter can be measured by averaging the ellipticities of an ensemble of 
galaxies. 



4.1. Distortion Matrix 



A distortion is an intensity-conserving mapping between the true position x and the observed 
position x' of the form 

x' = x'(x) = x + 5x(x). (16) 
The observed intensity i'(x) is thus related to the true intensity i(x) by 

iV) = *(*(*'))• (17) 
The local properties of the distortion are quantified by the distortion matrix <fi which is defined as 



« + 7i 72 + p \ _ d(5xj) 
72 - P « - 71 / 9xj 



(18) 
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where k is the convergence, 7$ is the shear, and p the rotation parameter. The convergence 
k describes overall rescalings, and the shear components 71 and 72 describe contractions and 
dilations parallel to, and at 45° from, the coordinate axes. In the case of weak lensing, <f> is directly 
related to the second-order derivative of the gravitational potential projected along the line of 
sight. While the rotation parameter p is expected to be negligible for weak lensing, it is generally 
not so for instrumental distortions. 

Let us consider a true position x = x° + Ax, which is close to a reference position x° (eg. 
the centroid of a source). For small distortions (i.e. if the offset Ax is small compared to the 
scale on which the distortion mapping varies), the corresponding distorted position is given by 
x\ = xf + {8ij + 4>ij) Axj + 0((/) 2 ), where xf = x^(x°) and 0{4> 2 ) denotes higher-order derivatives 
of the distortion mapping. Similarly, the true position corresponding to a distorted position 
x' = x 0/ + Ax', is given by Xj = x\ + (8ij — 4>ij) Ax' + 0(<p 2 ). By inserting this expression in 



equation (17) and by letting xf = x° = and Oi = Ax^, we obtain 

i'{e i )=i(e i -<i> ij e j ) + o{<t) 2 ). (19) 



4.2. Effect of Distortion on the Source Moments 



The quadrupole moment Z-j (Eq. Q) for the distorted image i'{9) (Eq. [19|) is given by 



4 = / d 2 e ei6jw{0)i(9 k - AiOi) + o{4> 2 ). (20) 



After Taylor expanding and integrating by parts we get, for a Gaussian weight function (Eq. B), 
I'ij = Iij + 21&[i0j-]fc + Iij4>kk — u~ 2 Iijki(t>ki + O(cj) 2 ), where the unprimed moments correspond to 
the undistorted moments. Similarly, the distorted monopole moment is related to the undistorted 
moments by I' = I + I(/)kk ~ ^^Ikifiki + 0(cfi 2 ). As a result, the distorted normalized moments 
(Eq. pi) are given by 

where the distortion susceptibility tensor is given by 



j; 7 = + D ijkl <j> kl + 0{4> 2 ) (21) 



Dijki = Dijki{3) = 25 k [iJj}i + u 2 (JijJki — Jijki) ■ (22) 

The brackets denote the symmetrizer, which, for an arbitrary tensor -A^ia—in °f ran k n i is defined 
by j 

A [hia-i n } = [ A hi2-i n + a11 n! permutations of {ii,i 2 , • • • (23) 



Inverting equation (21) yields 

J i j = J' ij -D' ijkl( j> kl + 0{ ( j> 2 ), (24) 

where D[j k i = Dij k i(J'). This expression can be used to correct the normalized moments for a 
known distortion. 
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4.3. Correction for the Fourth-Order Moments 

In principle, the correction for the fourth-order moments Jijki can be derived in a similar way. 
However, when the weight function w(6) is taken into account, the resulting expressions contain 
sixth-order moments and are very cumbersome. As we will see in §||, it is sufficient, in practice, 
to consider the unweighted (u> — * oo) corrections to Jijki- In this approximation, the corrected 
fourth-order moments Jijki are related to the distorted moments JL^l by 

Jijki - 4jkl ~ 4J m[ijfc^]m + 0((f> 2 ), (25) 

where, as before, the brackets denote the symmetrizer (Eq. p3fl ). This expression can be used to 
correct «/Ljy for instrumental distortions. 



4.4. Measurement of the Shear 

We now show how the gravitational shear can be measured by averaging over galaxy 
ellipticities. For this purpose, we consider galaxy moments that have been corrected for all 
instrumental effects, i.e. from instrumental distortion and PSF convolution, using the prescriptions 
presented in the other sections of this paper. The effect of a weak general distortion faj 
(Eq. |Igj]) on the (corrected) ellipticity of a galaxy can be derived by substituting the distorted 
moments J-j (Eq. ^l]) into the definition of the ellipticity e ? - (Eq. £L2|). This results in a relation 
between the distorted ellipticity e, t and the distortion parameters S{ = {k, 71,72; p} of the form 
€i = AijSj + 0(4> 2 ), with 

A _ ( Ml + Aei 2 - 2ef - t\H\ - \\ - \v\ -2eie 2 - e x /i 2 - 5^2 2e 2 \ , , 

ij ~ \ 112 + Ae 2 -2eie 2 - e 2 fn -\v 2 2 - 2e| - e 2 fi 2 - + \v x -2e x )' 



where A, \i% and U{ are the spinors defined in §3.2. Note that p is kept here for completeness, but 



is expected to vanish for weak lensing distortions. 

To measure the shear 73, we then average over an ensemble of galaxies which are assumed 
to be randomly oriented, intrinsically. Thus, all we need is the rotational average of the above 
relation, which we write as (e^) = {Ajj)5j + 0{<p 2 ). After discarding all terms which are not 
rotationally invariant, we obtain 



where 

Gi^2-(e 2 )-i(A)-i(e-/i), G^^x/,), (28) 

and e 2 = ef + e|, e • fi = t\p\ + £2/^2 > and e x \i = eifi 2 — £2/^1- Note that, to this order, the 
convergence k (and the rotation parameter p) do not affect the mean ellipticity (e t -). 
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Since we do not expect the galaxy population to have a preferred handedness, we can also 
discard terms which are not parity invariant. It is easy to check that e x [i is such a term since 
it changes sign when it is transformed to a left-handed coordinate system (eg. x — > x,y — ► —y). 
We are therefore left with a remarkably simple relation between the mean ellipticity (ej) and the 
shear 7, given by 

(et) = G 7i + O(0 2 ), (29) 
where G = Gi is the shear susceptibility. This expression agrees with equation (B13) in KSB, as 
corrected by HFKS. The inverse relation, 

li = E(e i ) + 0(c/ ) 2 ), (30) 
where E = G _1 can be used to measure the shear from the averaged ellipticity. 



4.5. Special Cases 

Unweighted moments: In this case, the effect of a general distortion fyj (Eq. |18|| ) on the 
ellipticity is more tractable. We find that the distorted ellipticity e[ is given by 

e • = ei + 2(<% - eiej)7j + 2e ij e jP + 0((j) 2 ), (31) 
where is the undistorted ellipticity, and the Levi-Civita symbol is defined by 
en = e22 = 0, ei2 = — e2i = 1. When averaged over an ensemble of randomly-distributed 
sources, this reduces to 

(^) = (2-(e 2 )) 7 , + O(0 2 ), (32) 
in agreement with equation (|29|) in the uj — > 00 limit. 



Circular source, unweighted moments: If we make the further simplifying assumption that 
the undistorted source is circular (e = 0), the distorted ellipticity (Eq. [31]) becomes 

e>^4 = 2 7l + O(0 2 ) (33) 

It is useful to plot the distortion ellipticity (rather than 7$) as a function of position on 
the chip, as a measure of the effect of the camera distortion (see figure |3[). 

Radial displacement: Let us consider the case where the displacement field (Eq. [p]]) is 
radial, i.e. where 5x = /(a;)§, where x = |x|, and /(x) is an arbitrary function. It is easy to 
show that, in this case, the distortion tensor (Eq. [Q) is 

, f df f\ XiXj f 

= U " x J — + (34) 
The corresponding distortion ellipticity (Eq. J33||) is 

df 



(35) 

where e x = {x\ — x|, 2#iX2}/(xf +x|) is the unit radial ellipticity field. From this expression, 
it is easy to see that will be radial (tangential) if — is positive (negative). 
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5. Convolution 

In this section, we study the effect of convolution by a weakly-anisotropic PSF on the source 
moments. The PSF can be decomposed into an isotropic and an anisotropic part. We show how 
each part can be corrected for, and how the PSF moments can be derived from stellar images. 



5.1. Effect of Convolution on the Source Moments 

Let us consider the case where the true galaxy image i(9) is convolved by a kernel p(9). We 
take p{9) to be normalized so that J d 2 6p(6) = 1, and centered so that / d 2 9 9ip{9) = 0. The 
observed image is given by 

if (e) = f d 2 e'p{e -e')i(0'), (36) 



so that the observed moments are 

= I d 2 e I d 2 e' e i e j w{e)p{e - e')i(o'). (37) 



lj — J "> " J "> » "t"j 

There are three angular scales in this equation, namely u, g, and a, corresponding to the 
size of the weight function w, the PSF kernel p, and the (unconvolved) source i, respectively. To 
simplify this expression, we need to make an expansion with respect to the ratio between two of 
these angular scales. We choose to expand with respect to ^, thereby assuming that the window 
function width is much larger than that of the PSF. In practice, the weight function scale is always 
chosen to be at least as large as the source size, i.e. co ^ a. To be conservative, we thus take a ~ to 
to collect the terms in the expansion. KSB instead, effectively expanded in powers of , 9 It 

V a +9 

is interesting to note that the expression we derive below for 1^ is nevertheless identical to theirs, 
to second order in — . 



After a change of variable {9" = 9 — 9'), and a Taylor expansion of w{9' + 9") about 9', the 
previous equation becomes I-j = hj + lP i3 - - ^I k ^P^ k - ^phjPkk + ^thjkiPki + (jjz) > where 
the brackets stand for the symmetrizer (Eq. ]23]|), and / and I^i are the (undistorted) weighted 
moments (Eq. As before, we take the weight function w(9) to be a normalized gaussian 
(Eq. |3|). The moments Pij are the (unweighted) moments of the convolution kernel, i.e. 

p.. = J d 2 9 9 i 9 ] p{9). (38) 

As stated above, this agrees, to this order, with KSB, as corrected for a factor of | by HFKS. 
In § |5.6| , we will show how the unweighted PSF moments Py can be derived from weighted 
stellar moments. In a similar fashion, it is easy to show that the convolved monopole moment is 
I' = I — -^jzIPkk + 2^IkiPki + o(^j^j. As a result, the convolved normalized moments (Eq. ||) 
are given by 

4- = Jij + C ijk iPki + O ( 4 J , (39) 
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where the convolution susceptibility tensor Cijki is given by 

2 1 

Cijkl = Cijkl(J) = Sik5ji — ^j2Jk[i$j]l + 7^4 [Jijkl — JijJkl] ■ (40) 

Inverting these equations yields 

ji3=4j-ci 3M Pk+o(^y (41) 

where C[^ kl = Cijki(J')- The last two equations can be used to correct the observed moments 
for the PSF, and only require knowledge of the observed quadrupole and fourth-order moments, 
Jij and Jijkl, of the galaxy, and of the quadrupole moments Pij of the PSF. We now show how 
this approximation can be used to correct for a weakly anisotropic PSF. 



5.2. Correction for the PSF 

In practice, it is often required to include galaxies with size a (and therefore weight-function 
width to) only marginally larger than the PSF size g. The PSF correction given by equation (S3), 
is thus not directly applicable, as the expansion series do not converge sufficiently fast. We can 
nevertheless apply the above correction scheme when the PSF is weakly anisotropic and sufficiently 
compact. In this case, we write the unweighted PSF moments as 

x-it)< (42) 

where g is the PSF radius, and is the PSF ellipticity, which is assumed to be small. We can 
then decompose the kernel p{6) into the convolution of an isotropic part p l {9) with an anisotropic 
part p a (9), as 

p = p i * p a . (43) 

It is easy to show that this implies that Pa = P/.- + P£, where P/,- and Ph are the unweighted 
moments of p l and p a , respectively. Without loss of generality, we further require that 

P\0) = ^e"^ (44) 

be a normalized circular Gaussian with standard deviation g. This implies that 

P?- = g%, and P§ = g 2 ^ | _ £ | J . (45) 

We now show how to correct each of these components in turn. 
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5.3. Anisotropic Correction 

If we consider a convolution with the anisotropic kernel p a (0) alone, the corrected moments 
(Eq. 0) become 

Jij = J'ii ~ C[ m I% + O (^Pj , (46) 

where, as before, C[j k[ = Cjjfei(J') is defined in equation (|40|). The residual terms are now 
supressed by the factor (e p ) 4 , and are thus negligible in practice. 



5.4. Isotropic Correction 

As explained earlier, the above approximation can not be applied to correct for the isotropic 
part of the PSF. However, since, by construction, the isotropic part p l of the kernel is a 
Gaussian, we can perform the isotropic correction exactly. Indeed, by inserting the form of p 1 
(Eq. Q|) into equation (j37|) and by integrating twice, we find the convolved moments to be 

/# = s£/ w «*y + (^) 4 iS 1 ' where 

ui^^+g 2 , g- 2 =g- 2 + u;- 2 , (47) 

and the superscript u and LUg in the moments I, Iij, and J^- indicate the standard deviation of their 
respective weight functions. The convolved monopole moment is simply I wl = I Ua . Consequently, 
the convolved normalized moments are 

$=(jj kj S'+9ZSv- (48) 

This convenient relation allows us to relate the convolved moments J?J with the unconvolved 
moments J^ 9 , but this time weighted with the wider weight standard deviation uj g . For later 
reference, we explicitly write the inverse relation 

which can be used to correct for the isotropic part of the PSF. 



5.5. Correction for the Fourth-Order Moments 

As in the case of distortions (§ |4.3|) , the general expression for the correction of the fourth-order 
moments Jijti for convolutions contains sixth-order moments and is very cumbersome. Here 
again, it is usually sufficient to consider the unweighted (u — ► oo) corrections to Jijki- With this 
approximation, the corrected moments Jijki are related to the convolved moments J-a and JLj^ by 

Jijkl ^ J'ijU ~ P m ~ § p [ijJ'kl] + QP [ij P kl], (50) 



-14- 



where 

p m = J d 2 e BiBjB^pie) (51) 

is the (unweighted) fourth-order PSF moment, and the brackets denote the symmetrizer (Eq. p3[ ). 
This expression can be used to correct J[j k i for convolutions. In the next section, we will show 
how Pijki can be estimated from stellar moments. 



5.6. Measurement of the PSF with Stars 

Stars are point-like and therefore have an intensity profile given by f(0) = S6W(0), where 
5 is the flux. The intrinsic moments (Eq. 0) of a star are thus I* = Sw*(0) and I*j = J*j = 0. 
We allow for the possiblity that the weight function, w*(0), for the stars be different from 
that for galaxies, w(0). After convolution (see Eq. JO] and above), the moments become 
/*' = /* - J^I*P kk + 0(4) and i?„' = I* Pa - J^I*P Akk + C>(fJ) , where P m was defined in 



2j(r Skk + O(g) and /*' = I*Pij - ^rP ijkk + 
equation ([5l"l). As a result, the normalized moments become 



T*/ 



P, 



ijkk 




(52) 



It is also easy to show that the observed fourth-order moments are given by I*j k[ 
while their normalized version is 



T*/ 

J ijkl 



Pijkl + O | 



i*p ijkl +o(ff), 



(53) 



We can invert these equations to obtain 



T*/ 



1 

2^1 



T*/ T*/ 

J ij J kk 



^ijkk 



+ o 



(54) 



and 



ijkl 



Jijkl + o 



(55) 



With these expressions, the unweighted PSF moments Pij and Pijki can be derived from the 
observed stellar moments J*- and J*- kk - These can then be corrected for the camera distortion 
using equations ( p4|p5| ) with uo — > oo, and then used in equations (46,49,^) to correct the galaxy 
moments. For this purpose, Pij and Pijki need to be interpolated across the chip. In practice, this 
can be done by fitting a low order polynomial to each component separately. Figure |l^ shows the 
PSF ellipticities for the WFPC2 camera derived from globular cluster observations (Figures |and 
||), after correction and interpolation. 

If the star and galaxy weight functions are equal {lo = u;*), then one only needs to keep the 
first term in equation (^). However, more accurate measurements of the stellar shapes can be 
achieved by taking a narrower weight function. In addition, it is desirable to avoid recomputing 



-15- 



the stellar moments for each value of u, which is often taken to be a function of the galaxy size. 
One would then choose uj* < lu, making equation (54) converge more slowly than equation (p]). 
In this case, one would thus need to keep the second term in equation ( |54| ) and to ensure that the 
residual error (o(jjf\) is acceptable compared to that for equation (p]) {oi^p. 



5.7. Special Cases 

Unweighted moments: In this case, the convolved moments become (Eq. [^] with u> — > oo) 
J-j = Jij + Pij. Thus, the observed square radius is 

d' 2 = d 2 + g 2 , (56) 

d 2 ti+g 2 



and the observed ellipticity is e- = dt+g** > where g 2 and ef are defined by equation ( |42| ) 

r deconv — j n 2s 

ij — J ij Si °ij 



It is sometimes useful to consider the moments J^ econv = J.. — g 2 5ij, which have been 



deconvolved from the isotropic part of the PSF only. The associated radius is d deconv = d, 

e deconv = e . + ( £ ) " 



while the associated ellipticity is e^ econv = a + (^) 2 ef. 



Unweighted moments, circular source: In this case, e% = and thus the observed ellipticity 
becomes 



while the isotropically deconvolved ellipticity becomes 

^deconv rdeconv J> I 9 \ p fzo\ 

e i = h e i - I 2 ) e i ■ i 58 j 

We have introduced the pre- and post-deconvolution reduction factors fd and f^ econv . These 
expressions are useful to estimate the effect of the PSF anisotropy on a galaxy. 



6. Combined Effect of Distortion and Convolution 

In practice, an image is deformed by a series of distortions and convolutions. The full 
treatment of the combined effect of a distortion and convolution is impractical given the 
complexity of the Dij^i and Cijki tensors (Eqs. [4^,22| ). However, the general behavior of the 



ellipticity is captured by considering the simplified case of a circular unweighted source which is 
subject to a weak distortion. This provides a good model of ellipticities observed in practice (see 
1| below). 

In the unweighted case (w(6) = 1) , the normalized quadrupole moments JL of an image 



deformed by a convolution followed by a distortion are (see Eqs. ]4l|24] with u> — > oo 



J'ij = {Jij + p ij) + {Jik + Pik)<t>kj + (Jjk + Pjk)<P ki + O(0 2 ), (59) 
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where Pjj are the PSF moments (Eq. [pq| ), ^ is the distortion matrix (Eq. [fig] ), and Jij are the 
undeformed moments. We will assume that the PSF is weakly anisotropic, i.e. that defined in 
equation ( |42| ) is small. Then, for an intrinsically circular source (Jy = d 2 5ij), the ellipticity of 
the deformed source is 



d 2 2 7t +g 2 (^ + 2 7 i 
ci 2 + 5 2 



+ O((0,e 



(60) 



where the deformation matrix faj was parametrized as in equation (|18|), and 5 is the PSF radius 
defined in equation (^). Notice that the convergence k and rotation parameter p do not appear 
in this first order expression. The observed square radius d' 2 = (J{ x + J^)/^ is given by 



d' 2 = (l + 2K)(d 2 + 5 2 )+O((0,ef) 2 ). 
For stars (d = 0), the observed ellipticity e* becomes 

e: = 6f + 2 7 , + O((0,6P) 2 ). 
We can thus rewrite e- in terms of observables as 



e' 



1- ( ^ 



27.+ (J) < + O((0,e 



(61) 



(62) 



(63) 



This expression is useful for computing the expected deformation of objects as a function of their 
observed size. This simplified model agrees well with the observed deformations in the Survey 
Strip (see §p|). 



7. Validity of the Method 
7.1. Shortcomings 

As was recently pointed out (Kaiser 1999; Kuijken 1999), the KSB method has several 
shortcomings^ which limit its validity and accuracy. First, the method requires a decomposition of 
the PSF into the convolution of an isotropic part and a compact anisotropic part. For most PSFs 
encountered in practice, this decomposition is formally ill-defined and the anisotropic part is not 
necessarily sufficiently compact. Another problem results from the fact that most PSFs do not 
fall off fast enough for the second moments (and higher moments) to converge. Consequently, the 
weighted second moments of the PSF depend strongly on the size of the window function. These 
problems are particularly severe for the PSF of HST which has broad wings extending beyond a 
central core, and thus cast doubts on the validity of the numerous weakdensing analyses based on 
HST observations. 



2 We thank Nick Kaiser, the referee, for pointing out and clarifying these shortcomings. 
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Since our method is based on the same principles as that of the KSB method, it also suffers 
from the same shortcomings. In addition, we have made the further assumption that, in the 
decomposition of the PSF, the isotropic part can be taken to be a gaussian (see § |5.2| ). This has the 
advantage of greatly simplifying the deconvolution, but can arguably produce further inaccuracies. 

These problems are certainly worrisome, and might eventually be solved by considering 
different methods such as those proposed by Kaiser (1999) and Kuijken (1999). However, the KSB 
approach has the advantage of being relatively simple and of requiring only a small amount of 
information about the PSF and galaxy shapes, namely the first few multipole moments. Moreover, 
these shortcomings are mainly formal in nature. In practice, one indeed always measures moments 
with a weight function width that is close to the width of the central core of the PSF, and 
which therefore do not diverge. For instance, for the HST PSF, the secondary sidelobes which 
compose the extended wings are smaller in amplitude than the central core by about one order of 
magnitude. Since the central core dominates in the convolution, one therefore expects that the 
corrections will be approximately correct. In the following paragraphs, we describe numerical 
simulations designed to test this assertion quantitatively. 

7.2. Simulations 

To test the weak link in our method, namely the PSF correction, we performed a series of 
numerical simulations. We first generated a WFPC2 PSF using the Tiny Tim software (Krist & 
Hook 1997). We chose the PSF to be particularly non-circular by placing it at the lower-right-hand 
corner of chip 2, at pixel coordinate (100,100) (see §|| for a description of the WFPC2 camera and 
of its PSF). We oversampled the PSF by using pixels which are 10 times smaller than the WFPC2 
pixels, and did not add any noise to the images. This isolates the systematic deconvolution errors 
from the random errors produced by pixelization and noise. To maintain uniformity, we will 
nevertheless quote angular sizes in WFPC2 pixels. 

We then measured the weighted moments Pij of the PSF as described in § |5.6| , using a range 
of weight function widths w*. For = 2 pixels, we found the PSF size and ellipticity components 
(Eq. |42|| ) to be g ~ 0.86 pixels, and ej ~ {—.009, —.048}. As expected, we found both g and e\ to 
diverge as increases. Fo the rest of the simulations, we set = 2 pixels, which is the stellar 
weight width that we used in our analysis of the survey strip (see § |8.3| ). 

We then convolved the PSF with elliptical gaussian "galaxies" of various intrinsic sizes d and 
ellipticities ej. Since the source of the problem lies with the PSF shape, this simplified galaxy 
model is sufficient for our purposes. We measured the moments of the convolved gaussian using 
a weight function width of uj = max(2, d) pixels, just as we did for the analysis of the survey 
strip (§ |3.3D . We then applied our PSF correction method (§[|) to obtain the corrected ellipticity 
gCorrected^ we ighted with an effective width uj g (Eq. [jj7| ). For an elliptical gaussian source, it is 
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easy to show that the true weighted ellipticity is e} Tue = ej 1 + (jj^) (1 — e 2 
induced by the PSF correction on the galaxy ellipticity is thus 



-i 



corrected 



Irue 



. The error Ae,- 



(64) 



We computed this residual error for a range of intrinsic galaxy sizes and ellipticities. The results 
of these simulations are presented in the next section. 



7.3. Results 

For the range of galaxy shapes which we considered (1.0 < d < 4.0 pixels and ej < 1), we 
found that, for a given galaxy size d, the ellipticity error Ae^ is very well approximated by 

Ae*~ Aei(0) + /(e)S, (65) 

where Ae«(0) is a constant, £i = e%/e is the unit radial ellipticity vector, and /(e) is a function 
describing the radial behavior. The constant term Ae«(0) is shown in figure [2|a for several relevant 
galaxy sizes, and has a modulus less than .004 for d > 1.5 pixels. Figure ^b shows the radial term 
/(e) as a function of the unweighted galaxy ellipticity modulus e, for the same galaxy sizes. The 
radial term depends strongly on the galaxy size d and ellipticity e, and has an amplitude of several 
percent for moderate ellipticities (e < 0.6) and sizes (d > 1.5 pixels). This amplitude is the size of 
the systematic error made in measuring the ellipticity of a single galaxy. 

However, for weak lensing measurements, one is not eventually interested in the ellipticity of 
a single galaxy, but in measuring a global change in ellipticity averaged over a galaxy ensemble. 
Let us consider a change in ellipticity ej — > e ■ = ej + 5{ produced by lensing^. The error in the 
measurement of dei is 

A^ = <Aei(0) * (Ae,) + 5e u (66) 

where the brackets refer to the average over the galaxy ensemble. For the functional form of 
equation (^5J), this becomes 

A6ei~ Aei(0) + h(e)5e h (67) 

where h(e) = \ /'(e) + . The first term was plotted on figure ^a. The second term is 
shown as a function of e and d in figure [|c, for an ellipticity change of 5e, = 0.05. Since, for our 
sample G ~ 1.4 (Eq. p9|), this value of 5e, corresponds to a shear of 7, ~ .04 and is therefore 



3 Strictly speaking, the ellipticity change Si produced by weak lensing has a weak dependence on a (See Eq.pl|). 
While this effect is only important for large ellipticities, the ensemble average is dominated by small ellipticities. For 
the purpose of estimating the size of the average correction error, we can thus ignore this small effect and consider a 
constant ellipticity shift. 
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representative of shear signals expected from large-scale structure. For moderate ellipticities 
(e < 0.6) and sizes (d > 1.5 pixels), the second term has an amplitude less than 0.004. 

Interestingly, both terms in equation ( |67|) tend to average out, when the distribution of 
galaxy sizes is considered. Indeed, in our treatment of the survey strip (see §|8.4|), we selected a 
sample of "large" galaxies with d > 1.5 pixels (and e < 1, which is not always satisfied because of 
noise). The average size for this sample is (d) — 2.58 pixels, which is close to the value of d for 
which both terms change sign (see figure ||). The ellipticity dispersion for this sample is a e ~ .31. 
Taking (d) and a e to be representative for this sample, we find the average residual errors to be 
A5ei ~ 0.001. This cancellation may, however, depend on the galaxy profile, while we have only 
considered gaussians in these simulations. To be conservative, we therefore take the residual error 
for our sample to be about 0.004, which is the error before averaging over d values. This is close to 
the residual error estimated from the anisotropic correction of the stars in our analysis of globular 
cluster fields (SQ. 

We also studied how the choice of the weight function widths uj and ll>* affect these results. 
We found that for moderate values (1 < < 3 and 0.5<i < uj < 2d), the magnitude of the residual 
errors did not change significantly from the above. However, the errors were much larger for uj 
and u;* outside of these ranges. This is, of course, a consequence of the divergences present in the 
PSF moments. 

We conclude that, the ellipticity error produced by the PSF correction can be several percent 
for an individual galaxy observed with WFPC2. However, it is only about 0.004 when averaged 
over a galaxy ensemble with d > 1.5 pixels, provided moderate values of the weight function 
widths are used. This does not alleviate the serious formal shortcomings of the method, and the 
need for a search for more robust methods (see Kaiser (1999); Kuijken (1999)). 

8. Application to HST images 

As an application of our method, we consider weak-lensing measurements with the WFPC2 
camera onboard HST. The camera consists of 3 800 x 800 pixel chips, with a pixel size of 0.1 arcsec. 
With the F814W filter, the PSF has a FWHM of about 0.09 arcsec. The instrumental effects for 
WFPC2 have been studied in detail by HFKS. They showed that the two main systematic effects 
for this instrument are the camera distortion and the PSF. We repeat their analysis by considering 
each of these effects in turn. We then apply these results to the galaxies found in the Survey Strip. 

8.1. Camera Distortion 

The camera distortion is caused mainly by the flat fielder, a refractive element in the WFPC2. 
This effect has been studied by Holtzman et al. (1995), who quantified this effect using astrometric 
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measurements of globular cluster observations. They modeled the distortion by a cubic polynomial 
in chip position for each chip. The polynomial coefficients can be used to derive a transformation 
matrix 09j imera between the WFPC2 chip coordinates (which are the distorted coordinates x' in 
Eq. [16 1) and the undistorted coordinates (x in the same equation). 



The distortion field is shown in figure ||. The quantity plotted is ef = 2ji, the ellipticity 
induced by the distortion field cpij on an intrinsically circular source (see Eq. |33]]). The pattern is 
radial and increases in magnitude with distance from the chip center. This differs from the results 
of HFKS who derived a tangential rather than radial camera distortion field. The radial nature 
of the pattern is confirmed by the fact that, using the coefficients of Holtzman et al. (f995) in 
equation (35), we find (J^ — > 0. (This can also be confirmed by drawing a box connecting 



the heads and tails of four adjacent arrows in Figure 15 in Holtzman et al.) 

To study the profile of the ellipticity pattern, it is useful to define the rotated ellipticity by 

ej = Rij(-2ip) ej, (68) 

where is the ellipticity of an object in the chip frame, Rij is the rotation matrix defined in 
equation [5], and (p is the polar angle of the source position about the center of the chip, measured 
counter-clockwise from the x-axis. A positive (negative) value of e\ corresponds to a radial 
(tangential) ellipticity pattern. A positive (negative) value of e£ corresponds to an anti-clockwise 
(clockwise) swirl pattern. The meaning of e\ is illustrated in figure ||. 

The profile of the rotated distortion ellipticity is shown as the dot-dashed line on figure 0. 
The mean ellipticity averaged over the 3 chips is listed in table |l]. The mean rotated ellipticities 
are found to be (e^) ^ 0.007 and (e£) < 0.001, while the mean absolute ellipticity (ej) is less than 
0.001. The effect of the camera distortion is thus small, but nevertheless comparable to the lensing 
signal expected from large-scale structure. 



8.2. Point-Spread Function 

The PSF is affected by diffraction by the telescope and scattering of light from different 
parts of the telescope and the optics. The PSF depends both on wavelength and time. The time 
dependence has two components (Krist and Hook, 1997). The first one is due to orbit-to-orbit 
"breathing" of the telescope. As the telescope orbits the earth every 90 minutes, it passes into and 
out of sunlight. The heating and cooling of the telescope changes its size, and thus its focus and 
PSF. The second time-dependent factor arises from a change in the focus of the telescope over 
longer periods of time. This change is produced by the outgassing of the graphite epoxy truss 
which supports the primary mirror. Approximately every six months, the secondary mirror is 
moved to bring the telescope back into optimal focus. The time-dependence of the focus position 
is summarized in Figure ||. 
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8.2.1. Tiny Tim 

As a first step, we used the Tiny Tim software to model the PSF of the HST (Krist & Hook 
1997). Tiny Tim takes into account the diffraction of light by the telescope to create a PSF, given 
the instrument (WFPC2 in our case), chip number, chip position, filter, spectrum of the object 
being modeled, and the focus of the telescope. Tiny Tim does not include the geometric distortion 
discussed above. 

We used Tiny Tim to model the PSF across the three WFPC2 chips. We chose the F814W 
filter and the spectrum of the objects to be similar to that of a type G star (B— V=0.619), which 
is typical of galaxies in the Survey Strip (see §|8.3j). We measured the moments (Eq. ||) and 
weighted them with a circular Gaussian (Eq. ||). For all stellar measurements, we adopted a 
weight-function width of uj* = 2 pixels, which resulted in an optimal sensitivity. 

To study the variation of the PSF across the chip, we used Tiny Tim to create a grid of PSFs. 
These PSFs were generated for a telescope with a mirror at a focus position of —1.5 microns, the 
value of the focus when the Survey Strip was taken (see §|8,3j). The resulting ellipticity pattern is 
shown in Figure |||. The corresponding ellipticity profile averaged over all chips is shown as the 
central dashed line in figure 0. In these figures, we have added the effects of the detector shear 
shown in figure [3|, so that a direct comparison can be made between these simulations and the 
globular cluster measurements which are presented below. The predicted PSF anisotropy is quite 
large ((e^) ~ 0.1 at the edge) and varies significantly across each chip and from chip to chip. 

To test the effect of focus changes, we created Tiny Tim PSFs at focus values of +2.0 microns 
and -5.0 microns, the full range of telescope focus (see figure |5|). The ellipticity profiles are shown 
are the upper and lower dashed curves in figure ^. The chip-averaged ellipticities for each focus 
value are listed in table |l|. Extreme focus changes from the central value of —1.5 microns produce 
variations of about 0.002 for (e'J) and about 0.006 for (ei) and (£2), when averaged over all chips. 

8.2.2. Globular Clusters 

For a more direct measurement of the PSF, we obtained archival images of the globular 
cluster M4 observed by Richer et al. (1995). We also obtained archival images of the globular 
cluster NGC6572. This was the method of calibration used by HFKS. It is fortuitous that the M4, 
NGC6572, and the Survey Strip data all have focus values of approximately —1.5 microns (see 
Figure |5|) . 

The M4 data consisted of 3 different fields in each cluster. Each field had 8 different exposures: 
four pairs of dithered pointings separated by very nearly 2 pixels. Other dithers were separated 
by non-integer numbers of pixels. We shifted the dithers which are separated by integer numbers 
of pixels and compared them pairwise to detect cosmic rays. A pixel was flagged as a cosmic ray 
if its flux was significantly above that of the corresponding pixel in its pair image. We selected 
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stars from the M4 exposures by excluding saturated stars and those that are faint and dominated 
by noise. We also excluded stars confused by close neighbors, and stars containing a cosmic ray 
flagged pixel within a 5 pixel radius. Our resulting sample consisted of about 160 stars per chip. 
We then averaged the three moments Pki measured in each of the eight exposures. The analysis of 
the NGC6572 data was simplified by the fact that we had only one field and the images were not 
dithered. 

The ellipticity fields for M4 and NGC6572 are shown in figures || and |9], respectively. The 
ellipticity profile for each globular cluster is shown as the upper and lower solid line in figure [7|, 
respectively. The chip-averaged ellipticities for each globular cluster data set are listed in table |]. 
The M4 and NGC6572 ellipticity patterns are qualitatively similar, and both show tangential 
ellipticities of the order of 0.1 at the edge of all chips. The M4 and NGC6572 average ellipticities 
(e^) are —0.040 and —0.053, respectively, and are both larger than that predicted by Tiny Tim 
(about —0.035). The fact that the mean tangential ellipticity for NGC6572 is larger than that for 
M4 by about 0.013, shows that the PSF varies substantially over periods as short as a few days. 
This change is larger than that predicted by Tiny Tim for extreme focus changes (see § |§.2.1 ). 



Moreover, the globular cluster ellipticities are qualitatively different from that predicted by Tiny 
Tim. For example, figure |l(] shows a "teardrop" shape to the M4 ellipticity pattern of chip 2, 
with a cusp at the top of the chip. In contrast, the Tiny Tim ellipticity pattern of the same chip, 
shown in figure || has a cusp on the right-hand side. Small changes in this cusp result by changing 
the focus of the telescope in Tiny Tim, but it does not change the qualitative pattern of the 
ellipticities. 

Because of these shortcomings, we used the observed globular cluster moments for the 
corrections, rather than the Tiny Tim predictions. To reduce patchiness, we combined the M4 
and NGC6572 stars, and corrected the moments for the weight function (Eqs. [ p4"| , |55| ] ) and for the 
camera distortion (Eqs [24,25| with u> — ► oo) to obtain the PSF moments P{j and Pijki- We then 



fitted a fifth order polynomial in chip position for each moment component and for each chip. 
Because the star counts in the two globular cluster are close to one another (i.e. about 160 per 
chip), this amounts to giving equal weight to each of them. The ellipticity profile for the combined 
globular clusters is shown as the central solid line in figure 0. Figure 10 shows the resulting fitted 
ellipticity pattern of the PSF. 

The mean PSF radius for the combined M4/NGC6572 was found to be g ~ 0.89 pixels (see 
Eq. fl42[). This is much smaller than the formal PSF radius that one would measure using a Tiny 



Tim PSF with CJ* — > oo. As explained in §[741, the quadrupole moments indeed formally diverge 
because of the extended wings of the HST PSF. Our value is however commensurate with the PSF 
FWHM (0.9 pixels), and should thus be considered as an effective PSF radius. 
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8.3. Survey Strip 

The Survey Strip is a set of 28 contiguous pointings with the HST in two colors, V (F606W) 
and I (F814W) (Groth et al. 1995; Rhodes 1999). The images were taken in March and April 
1994 with the WFPC2. The Strip has about 10,000 galaxies down to / ~ 26 and covers an area of 
about 108 square arcmin over a 3'. 5 x 44'. region. The Strip has already proved a useful data set 
in exploring the distribution of luminous matter through number counts and the 2-point angular 
correlation function (Rhodes, Groth and Refregier 1998; Rhodes 1999). We are in the process of 
exploring the large scale distribution of matter through weak lensing of the Strip galaxies (Rhodes 
1999; Rhodes, Refregier, & Groth 1999). 

Here, we use the Strip galaxies to test our shear measurement method. For this purpose, we 
used the I images only. The catalog of galaxies was created using the Faint Object Classification 
and Analysis System (FOCAS) within IRAF (Jarvis &; Tyson 1981). An object was considered 
as detected if two contiguous pixels were more than 3a above the sky background. However, an 
object was included in the final catalog only if both the I and V images had a detection within an 
error box of approximately 5 pixels (0.5 arcsec). This resulted in a sample of 9448 galaxies with 
I < 26. 

The moments of the galaxies were calculated as described above (Eq. || ) . The variable weight 
function was chosen to be u> = max(2, y^"), where A is the detection area calculated by FOCAS. 
The lower bound of u was chosen to match the width used for stars. 

The importance of the camera distortion and of the PSF depends on the radius d of the galaxy 
(Eq. pT|). We thus considered subsamples of small and large galaxies, with 1.0 < d! < 1.5 pixels 
and d! > 1.2 pixels respectively. Here d' is the root mean squared observed radius (Eq. HH) and 
g ~ 0.89 pixels is the PSF radius (see § 8.2. 2| ). Because of the noise, a fraction of the galaxies have 



unphysical ellipticity values (e > 1). Because these outliers would dominate ellipticity statistics, 
we only retained galaxies with e < 1. The mean observed radius and reduction factor (Eq. j57j) 
are (d 1 ) ~ 1.26 pixels and fd — 0.50, and {d') ~ 2.16 pixels and fd — 0.17, for the small and large 
galaxies, respectively. Each subsample respectively comprises 41% and 77% of the total number of 
/ < 26 galaxies. 



The ellipticity profile for each subsample is shown as a solid line in figures [11] and 12. Also 
shown are the asymptotic cases, namely the combined M4/NGC6572 stars (d = 0), and the 
camera distortion ellipticity (d — > oo). The dashed line shows the prediction of the simplified 



model described in §|6] (see Eq. [63]). The small galaxy sample shows a tangential ellipticity which 
increases with radius from the chip center, in agreement with the model prediction. The large 
galaxies, on the other hand, do not display any significant tangential ellipticity, again in agreement 
with the model. In both cases, (e^) is consistent with zero. 



Using the camera distortion matrix derived in §|8.1| and the PSF moments derived in §8.2, we 



then corrected the Strip galaxies. The ellipticity profiles at different stages of the correction for 



-24- 



the small and large galaxy subsamples are shown in figures and [[I]. Since the camera distortion 
is radial, its correction increases the tangential ellipticity of the galaxies. On the other hand, 
the correction for the PSF, which is mainly tangential, reduces the tangential component of the 
ellipticity. For both subsamples, the resulting corrected profile is consistent with 0. As expected, 
the PSF correction has a smaller effect on the large galaxies than on the small galaxies. 



8.4. Error Budget 

The effect of the corrections for the chip-averaged ellipticities are summarized in table [l]. To 
ensure statistical independence, we considered in this table a subsample of large galaxies with 
d' > 1.5 pixels (and e < 1). This subsample contains 51% of the total number of galaxies with 
I < 26, corresponding to a surface density of n — 32 arcmin -2 . Their mean observed radius 
is (d') ~ 2.58 pixels, corresponding to a post-deconvolution reduction factor of f^ econv ~ 0.13 
(Eq. @|). As a test of the correction algorithm, we also corrected the globular cluster stars, 
exactly as we corrected the galaxies. (We have not corrected either the stars or the small galaxies 
for the isotropic PSF, as this would produce diverging ellipticities). The results for the corrected 
stars are also shown in the table. 

The self-corrected M4/NGC6572 stars have a residual ellipticity of (ej) ~ .003 and 
(ei) ~ (£2) — .002. This is a measure of the errors resulting from the approximations in our 
correction method and from the fit to the stellar moments. Another estimation of the errors 
from the method was presented in §[?]. By performing numerical simulations, we showed that the 
residual ellipticity errors in the PSF correction was 0.004, after averaging over a galaxy ensemble 
with d > 1.5 pixels. Given the agreement of these two estimates, we take the error in the correction 
method to be about 0.004. 

As we noted in § |8.2j , the M4 and NGC6572 ellipticities differ by about 0.01 for both rotated 
and absolute ellipticities. After the anisotropic PSF correction, the small galaxies have a residual 
ellipticity of (e^) — 0.004, (£2) — -011 and (ei) ~ .005. Since the small galaxies have a reduction 
factor of fd — 0.50, their residual (ei) ellipticity indicates that the PSF variation is about 0.02, 
which is larger than that given by the comparison of M4/NGC6572. We therefore take the 
variations of the chip-average PSF ellipticity to be about 0.02. 

How does this uncertainty affect the ellipticities of the large galaxies? As noted above, the 
post-deconvolution reduction factor for the large galaxies (with d' > 1.5 pixels) is f^ econv ~ 0.13. 
Consequently, our uncertainty of 0.02 in the PSF ellipticity produces an uncertainty in the 
ellipticities of the large galaxies of only 0.003. This is of course a commendable consequence of 
the small size of the WFPC2 PSF. Note also that, for the large galaxies, the camera distortion 
correction is about the same size as that for the PSF, so that the two corrections almost cancel 
each other. The changes in ellipticities produced by each of these corrections is less than 0.005 for 
the large galaxies. The total systematic uncertainty is a combination of the correction uncertainty 
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and of the PSF variability, and is thus about 0.004. 

The statistical uncertainties are determined by the ellipticity variance c 2 = ^(e 2 ) = 
\ (( e i) + ( e i))- Both noise in the image and the intrinsic shapes of the galaxies contribute to this 
dispersion. The rms ellipticity a t is listed in the last column of table |l], at various stages of the 
correction. For both the large and small galaxies, a e does not increase when we correct for the 
camera distortion and for the anisotropic PSF. This shows that our correction method does not 
introduce any appreciable noise in the ellipticity measurements. Not surprisingly, a e increases 
moderately when we correct for the isotropic PSF. This is expected since this deconvolution 
reduces the galaxy size, and thus reduces the denominator in the definition of the ellipticity. For 
the large galaxies, the rms ellipticity after all corrections is about a e ~ 0.31. Since the mean 
number of such galaxies per chip is about N g ~ 57, the la sensitivity to detect the shear in a 
chip is a e /yf{Ng) ~ .04. This is close to the expected rms shear expected from weak lensing by 
large-scale structure in 1.3 arcmin cells, for cluster normalized CDM models and for a source 
redshift z = 1 ( Jain & Seljak 1997). This shows that the expected signal-to-noise ratio of a single 
WFPC2 chip is about 1. 



9. Conclusions 

We have revisited the KSB method to measure the weak lensing shear from the shapes of 
galaxies. In our method, the corrections for the camera distortion and PSF convolution are 
performed using moments rather than ellipticities. Using a gaussian weight function, we derived 
explicit expressions for the corrections, which involve only second and fourth order moments. 
We clarified the convergence of some of the approximations made by KSB, and showed how the 
weight function for stars can be chosen to be different from that for galaxies. We also showed how 
the isotropic part of the PSF can be assumed to be a gaussian, to the required level of precision, 
and can thus be corrected for exactly. We derived the explicit relation between the shear and the 
ellipticities by decomposing moments into tensors with definite rotational properties. 

We addressed the recently exposed shortcomings of the KSB method (Kaiser 1999; Kuijken 
1999). Our method, as well as the KSB method, has formal problems arising from the fact 
that PSFs encountered in practice are not sufficiently compact. We used numerical simulations 
to assess the importance of these problems in the analysis of WFPC2 images. We found that 
the ellipticity error produced by the PSF correction can be of several percent for an individual 
galaxy. However, it is only about 0.004 when averaged over a galaxy ensemble with d > 1.5 pixels, 
provided moderate weight function widths are used. 

We studied systematic effects arising in WFPC2 images. From globular cluster observations, 
we confirm the results of HFKS, who found ~10% PSF ellipticities at the edge of each chip. We 
find however that the camera distortion is radial, rather than tangential. It produces ellipticities 
of the order of 0.7%. We further find that the PSF ellipticity varies by as much as 2% over time. 
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We applied our correction method to the HST Survey Strip. We showed that the different 
stages of our correction do not introduce any appreciable noise. We studied the dependence of 
galaxy ellipticities on the galaxy size. Small galaxies are more sensitive to the PSF and also 
indicate that the PSF varies with time. For large galaxies (observed radii d! > 1.5 pixels), the 
total systematic uncertainty is about 0.4%, and results from a nearly equal contribution from the 
correction uncertainty and from the PSF variability. The statistical la uncertainty in measuring 
the shear in a single WFPC2 1.3 x 1.3 square arcmin chip is about 4%, for this subsample of 
galaxies. This provides good prospects for detecting a cosmic shear signal with the strip and other 
deep HST surveys. In Rhodes et al. (1999) and Rhodes (1999), we will describe our search for 
such a signal. 
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Table 1. Effect of Corrections on Chip- Averaged Ellipticities 
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Corrections for: d: camera distortion, a: PSF anisotropy, i: isotropic PSF. The PSF is derived from the combined 
M4+NGC6572 stars. 

V = i< £ 2 ) = i((e?> + (el)) 

C PSF predicted by Tiny Tim for several focus values 

d Galaxies with magnitudes I < 26, ellipticities e < 1, and observed radii 1.0 < d' < 1.5 and d! > 1.5 pixels, for 
the small and large sample, respectively. 
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Fig. 1.- 



Synopsis of the weak lensing measurement method. 
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Fig. 2. — Residual errors of the PSF correction determined from numerical simulations. Panels 
(a), (b) show, respectively, the constant term Aej(0) and the radial function /(e), for measuring 
absolute ellipticities. Panel (c) shows the radial function h(e) for measuring an ellipticity change of 
5ei = 0.05 (see text). These quantities are plotted as a function of the unweighted galaxy ellipticity 
e and size d (in WPC2 pixels). 
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Fig. 3. — The WFPC2 Detector Shear pattern. The quantity plotted is ef = 2ji which is the 
cllipticity induced by the detector shear for an intrinsically circular source. Chips 2, 3, and 4 are 
in the upper-left, lower-left, and lower-right corners, respectively. Chip 1 (the Planetary Camera) 
is not shown. Each chip is 1.3 arcmin on a side. 
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Fig. 4. — PSF ellipticity predicted by Tiny Tim. This figure corresponds to a focus of — 1.5/xm, close 
to the mean value for the Survey Strip observations. The ellipticities are shown after application 
of the camera distortion. 
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Fig. 5. — HST focus as a function of time. The vertical lines give the times of mirror movements 
generated to bring the telescope back into focus. The diagonal lines represent the average mirror 
position as a function of time. Orbit-to-orbit "breathing" produces considerable variations (up to 
several microns) about these mean values. The rms value of these variations is shown by the dotted 
lines. This figure is patterned after Figure 5.7 in Biretta (1996). The values corresponding to the 
observations of the Survey Strip, and of the M4 and NGC6572 globular clusters are also indicated. 
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Fig. 6. — Illustration of the meaning of rotated ellipticity e\. This ellipticity is defined by choosing 
coordinates axes which are rotated about the chip center. The patterns corresponding to positive 
and negative values for each component of e[ are shown. 
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Fig. 7. — Rotated ellipticity profile for the PSF and the camera distortion, as a function of 
radius about the chip center. The solid lines show the ellipticity as measured for M4, combined 
M4/NGC6572 , and NGC6572 stars, from top to bottom respectively. The dashed lines show the 
PSF ellipticity predicted by Tiny Tim for focus values of +2.0, —1.5, and — 5.0^m, from top to 
bottom respectively. The dot-dashed line shows the distortion ellipticity ef . 
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Fig. 8. — Ellipticities of the stars in the globular cluster M4. A Gaussian weight function of width 
uj* = 2 pixels was used to measure the stellar moments. 
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Fig. 9. — Ellipticities of the stars in the globular cluster NGC6572. 
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Fig. 10. — PSF ellipticities derived from the combined M4/NGC6572 stars after being corrected 
for shear and weighting and interpolated. 
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Fig. 11. — Rotated ellipticity profile for small galaxies in the survey strip. Galaxies were selected 
to have I < 26, and radii 1.0 < d! < 1.5 pixels. The PSF and the camera distortion ellipticity ef ,r is 
also shown for comparison. The profile for the Strip galaxies expected in the simplified theoretical 
model is also shown (see text). 
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Fig. 12. — Same as the previous figure, but for large galaxies (d! > 1.2 pixels). 
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Fig. 13. — Effect of the corrections for small galaxies (1.0 < d! < 1.5 pixels, / < 26) in the survey 
strip. The profile for the galaxies is shown at different stages of the correction algorithm. The error 
bars are similar for each of the three profiles, but, for clarity, were only displayed for the last case. 
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Fig. 14. — Same as the previous figure, but for large galaxies (d! > 1.2 pixels). 



